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1 Introduction: gluon saturation and structure functions 

r-| ! At very high energies corresponding to very small values of Bjorken x variable 

the density of partons in the hadronic and nuclear wave functions is believed 
to become very large reaching the saturation limit (see [T|2] and references 
therein). In the saturation regime the growth of partonic structure functions 
with energy slows down sufficiently to unitarize the total hadronic cross sec- 
tions. The gluonic fields in the saturated hadronic or nuclear wave function are 
very strong p!p] . A transition to the saturation region can be characterized 
by the saturation scale Ql{s), which is related to the typical two dimensional 
density of the partons' color charge in the infinite momentum frame of the 
hadronic or nuclear wave function. The saturation scale is an increas- 

ing function of energy s and of the atomic number of the nucleus A. At high 
enough energies or for sufficiently large nuclei the saturation scale becomes 
much larger than Aqqj^ allowing for perturbative description of the scattering 
process at hand. The presence of intrinsic large momentum scale Qs justifies 
the use of perturbative QCD expansion even for such a traditionally non- 
perturbative observables as total hadronic cross sections. There has been a 
lot of activity devoted to calculating hadronic and nuclear structure functions 
in the saturation regime. The original calculation of quark and gluon distri- 
bution functions including multiple rescatterings without QCD evolution in 
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a large nucleus was performed in [3]. The resulting Glauber- Mueller formula 
provided us with expressions for the partonic structure functions which reach 
saturation at small Q^. McLerran and Venugopalan has argued that the large 
density of gluons in the partonic wave functions at high energy allows one to 
approximate the gluon field of a large hadron or nucleus by a classical solu- 
tion of the Yang- Mills equations. The resulting gluonic structure functions has 
been shown to be equivalent to the Glauber-Mueller approach p!p] . 

To include quantum QCD evolution in this quasi-classical expression for the 
structure functions one has to resum the multiple BFKL pomeron exchanges. 
The evolution equation resumming leading logarithms of energy (a^ In s) and 
the multiple pomeron exchanges was written in |lj using the dipole model of [5] 
and independently in [6J using the effective high energy lagrangian approach. 
The equation was written for forward scattering amplitude A^(£, 1^) of a 
quark-antiquark dipole with transverse size r at impact parameter h with 
rapidity Y scattering on a target hadron or nucleus, which, in turn can yield us 
the F2 structure function of the target which is measured in DIS experiments. 

The evolution equation for closes only in the large-A^^c limit of QCD and 
provides the basis for all phenomenological applications. It describes the onset 
of the gluon saturation for small x and/or larger A. The most striking con- 
sequence of the gluon saturation is scaling of the total DIS cross section at 
small X with variable the so-called geometric scaling, observed 

in DIS on both proton [7] and nuclear targets||8]. 

Several other observables can be calculated in the framework of the saturation 
approach to hadronic and nuclear collisions. In this paper we review applica- 
tions of the saturation approach for calculation of various inclusive observables 
in pA collisions. 



2 Inclusive gluon production 



Single inclusive gluon production in the quasi-classical approximation has been 
derived in [H]. The inclusive cross section for the scattering of a dipole of 

transverse Xqi on the target reads 



TT J i,j=0 1—1 — «l 1—2 — il 



d'^k dy 
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This formula sums up the multiple rescatterings of the qq pair and the pro- 
duced gluon on the nucleons in the target nucleus in the = light cone 
gauge (+ is the direction of motion of the incident color dipole). Inclusion of 
quantum evolution amounts to resuming all possible real and virtual gluon 
emissions in addition to the emission of the measured gluon. In spite of the 
explicit breaking of factorization for individual diagrams in any known gauge 
the sum over all diagrams can be written in a simple fcr-factorized form [TO] 

«*=c;f/''*'<2M^<^-2). (2) 

where the unintegrated gluon distribution function is defined as 

Hx,k') = / dVre-'^-^-VlNG{r,b,y = ln(l/x)) , (3) 



with NG{r,b,y) the forward amplitude of a gluon dipole of transverse size r 
at impact parameter b and rapidity y scattering on a nucleus. If the forward 
scattering amplitude of qq dipole is found from the evolution equation, then 
the gluon dipole scattering amplitude can be calculated from 

NG{r,b,y) = 2Nir,b,y)-N'ir_,b,y). (4) 



Let us emphasize that in this approach we consider gluon saturation effects 
only in nucleus, while treating proton as a dilute object. Analysis of the in- 
clusive gluon production has been also performed in [TT|12fl3j . 

This approach to the high energy pA collisions leads to a rather successful 
phenomenological applications which allow the direct comparison with the 
experimental data, see Fig. [TJ While the multiple rescatterings in the nucleus 
without quantum evolution lead to mere redistribution of the gluon's trans- 
verse momentum towards higher values (Cronin effect), the effect of quantum 
evolution is to tame the growth of the scattering amplitude at higher ener- 
gies/rapidities and/or for heavy nuclei. As the result the saturation approach 
predicted the onset of suppression of the nuclear modification factor in dAu 
collisions at forward rapidities P^15lll6fl3fl7j . see Fig. [H 

It was suggested already in the pioneering papers on the gluon saturation 
that the particle correlations in the saturation regime must be significantly 
suppressed as compared to the low density regime (see [2]). Indeed, unlike 
in hard processes where two jets are produced back-to-back, in saturation 
regime most gluons are produced with a semi-hard momentum of the order 
of Qs ( "monojets"). From the McLerran-Venugopalan model point of view, 
the emitted gluons correspond to the commuting classical non-Abelian fields. 



3 




(a) (b) 

Fig. 1. (a) Hadron multiplicity in d-Au collisions |18|19j , (b) nuclear modification 
factor as a function of the transverse momentum px at pseudo-rapidity ry = 3.2[20]. 

The double inclusive gluon production in the quasi-classical approximation has 
been addressed in [21 22] where it was argued that the monojet correlations are 
responsible for generation of the azimuthal asymmetry in heavy ion collisions 
at large pr- Quantum evolution effects where discussed in [23] where it was 
predicted that suppression of correlations at forward rapidities is larger than 
at central ones and that this suppression grows with the rapidity interval 
between the produced hadrons, see Fig. [21 Other possible signatures of the 
gluon saturation in the hadron correlation function were discussed in j24j . 
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Fig. 2. (a) Azimuthal correlations between forward and backward hadrons, (b) scaled 
elliptic flow variable f 2 • M as a function of px [25] . M is multiplicity. 
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Exact theoretical result for the double-gluon production case has been ob- 
tained in [26]. It turns out that the /cy-factorization fails. Instead a more 
complicated factorization picture emerges. 



3 Heavy quark production 

Heavy quark production in hadronic collisions in high energy QCD is one of 
the most interesting and difficult problems. It is characterized by two hard 
scales: heavy quark mass m and the saturation scale Qg- The threshold for 
the invariant mass of the quark q and antiquark q production is 2m. There- 
fore, if m is much larger than the confinement scale Aqcd, it guarantees that 
a non-perturbative long distance physics has little impact on the quark pro- 
duction making perturbative calculations possible (for a review see [27]). For 
all processes involving heavy quarks with momentum transfer of the order of 
~ large saturation scale implies breakdown of the coUinear factorization 
approach. The factorization approach may be extended by allowing the incom- 
ing partons to be off-mass- shell. This results in conjectured fcr-factorization. 
Although the phenomenological applications of the fcr-factorization approach 
seem to be numerically reasonable at not very high energies [28] its theoreti- 
cal status is not completely justified. Like coUinear factorization it is based on 
the leading twist approximation. However, at sufficiently high energies, higher 
twist contributions proportional to {Qs/mY^ become important in the kine- 
matic region of small quark's transverse momentum, indicating a breakdown 
of factorization approaches. 

The fact that the saturation scale at high enough energies and for large nu- 
clei is large, Qs ^ ^qcd, combined with the observation that the typical 
transverse momentum of particles produced in pA scattering is of the order 
of that saturation scale, leads to the conclusion that Qg sets the scale for 
the coupling constant, making it small. This allows one to perform calcula- 
tions for, say, gluon production cross section in pA collisions using the small 
coupling approach [lOj. The same line of reasoning can be applied to heavy 
quark production considered here: the saturation scale Qg is the important 
hard scale making the coupling weak even if the quark mass m was small. 
Having the quark mass m as another large momentum scale in the problem 
only strengthens the case for applicability of perturbative approach. 

Production of quark-antiquark pairs in high energy proton-nucleus collisions 
and in DIS has been calculated in [2^30|31II32II33] . The results of the calcu- 
lations in a model based on the /cT-factorization approach are displayed in 
Fig. [3](a). The heavy quark mass delays the onset of the gluon saturation ef- 
fects to higher rapidities ~ 2 (for charm). The recent experimental results 
for the nuclear modification factor for muons associated with D-meson decays 
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are consistent with this prediction, see Fig. [3](b). 
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Fig. 3. (a) Charmed meson yield as a function of iVcou 
factor for muons associated with D-meson decays 
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4 Future perspectives 



In this admittedly brief review article we discussed how the effect of gluon 
saturation can be taken into account in the framework of the perturbation 
theory in the case of gluon and quark spectra, multiplicities and correlations. 
These observables are in agreement with the recent experimental data. In ad- 
dition to the aforementioned calculations of the gluon and quark spectra and 
multiplicities an extensive work has been done on computing the production of 
J/\Ef[36], valence quarks |4T]|42] . prompt photons[37] and di-leptons [38|39f40j . 
These observables still require detailed experimental investigation. Although 
we observed the onset of the gluon saturation at RHIC, this effect will be far 
more dramatic at LHC In fact, production of bulk of particles in AA and pA 
collisions at midrapidity at LHC is predicted to be governed by the satura- 
tion regime. Measurement of the energy dependence of the various physical 
quantities is an important test for the gluon saturation approach. 
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